Quadrupolar correlations and deformation effect on 
two neutrino £j3^ and ee modes of ^^^Dy isotope 



P K Rathi, R Chandra! 2, g singh\ P K Raina^ and J G 
Hirsch^ 

^Department of Physics, University of Lucknow, Lucknow-226007, India 
^Department of Physics and Meteorology, IIT, Kharagpur- 721302, India 
^Institute de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
A.P. 70-543, Mexico 04510 D.F., Mexico 
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are studied in the Projected Hartree-Fock-Bogoliubov framework for the 0"*" — > 0"*" 
transition. Theoretically calculated half-lives of electron-positron conversion and 
double-electron capture modes are presented. The effect of the quadrupolar 
deformation on nuclear transition matrix element M2iy is also investigated. 
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1. Introduction 

The nuclear decay, which is mediated by strangeness conserving charged weak 
current is expected to occur through double-electron emission {(3~(3~), double-positron 
emission (/3+/3+), electron-positron conversion and double-electron capture (ee) 

with the emission of two neutrinos, no neutrinos, single Majoron and double Majorons. 
The P'^(3~^ /e(3^ /ee processes are energetically competing and we refer to them as e+/3/9 
decay. The (/3/3)2^ decay conserves the lepton number exactly and is an allowed process 
in standard model of electroweak interactions (SM). The lepton number violating 
{PP)q^ decay, which is theoretically possible in many a gauge models beyond the SM, 
has not been observed so far. All the present experimental efforts are devoted to 
its observation, which would immediately imply that neutrinos are massive Majorana 
particles. In comparison to the P~ fi~ decay, less attention has been paid to study the 
e'^ 13(3 decay due to relatively low Q-values and low abudances of e^ (3(3 emitters. The 
experimental and theoretical developments in the study of nuclear (3(3 decay have 
been excellently reviewed over the past years [Tl[l[l[l[Sl[ni[71[Hl[ni[TU]. 

The half-life T'^j2 °^ {PP)2u decay is a product of accurately known phase space 
factor G21/ and nuclear transition matrix element (NTME) Hence, the validity 

of different models employed for nuclear structure calculations can be tested by 
comparing the calculated NTMEs M2V with those extracted from the experimentally 
observed half-ives Tfj2 ^'^^ ~* transition of (/3~/3~)2^ mode. In contrast 

to the {(3~(3~)2^ mode, which has been observed in ten (3~ (3~ emitters out of 35 
potential candidates, the (e+/3/?)2j, decay modes have not observed experimentally so 
far. However, limits on their half-lives T'^j2 have been given for 24 out of 34 possible 
isotopes ^llj. Meshik et al have reported on a positive observation of ^^°Ba decay with 
half-life T1/2 = (2.16±0.52) x 10^"'^ y for all modes [T^ and this result is consistent with 
theoretical expectations for the (e£)2^ mode TS'. This value is in slight contradiction 
with the experimental limit T1/2 > 4 x 10^^ y reported in [M]. If confirmed, it would 
be the very first observation of e+/3/3 decay. In the absence of experimental data, 
there is no way to judge the reliability of present nuclear structure calculations. In 
principle, the (3~ P~ decay and e+/3/3 decay can provide the same information. The 
observation of (e+/3/3)2j^ decay will be interesting from the nuclear structure point of 
view, as it will be a challenging task to calculate the NTMEs of these modes along with 
the {(3~ (3~)2jj mode in the same theoretical framework and it will further constrain 
the nuclear models employed to study the (3(3 decay. In addition, the observation 
of (e^/3/9)Q^ decay will be helpful in deciding finer issues like dominance of mass 
mechanism or admixture of the right handed current in the electroweak interaction 

Rosen and Primakoff were the first to study the {e~^(3f3)2^ decay theoretically 
[1]. Later on, Kim and Kubodera estimated the half- lives of all the three modes with 
estimated NTME and non-relativistic phase space factors [15]. Abad et al performed 
similar calculations using relativistic Coulomb wave functions [16| . In addition, some 
other theoretical studies were done by Zel'dovich et al [17], Eramzhyan et al [18], 
Bernabeu et al [19 and Balaev et al [20]. The {e'^ (3(3)2^ decay has been studied 
mainly in three types of models, namely shell model and its variants, the quasiparticle 
random phase approximation (QRPA) and its extensions and alternative models |7]. 
The nuclear many-body problem is solved as exactly as possible in the shell-model 
and it is the best choice for the calculation of the NTMEs. However, the number of 
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basis states increases quite drastically for most of the Pf3 decay emitters as they are 
medium or heavy mass nuclei and hence, Vergados has studied the (e+/3/9)2^ decay of 
^^Ni, ^^Mo and ^^Ru nuclei in the weak coupling limit. Over the past years, the large 
scale shell-model calculations have been successfully performed to study the potential 
PP emitters [2TJ [H [23l [24] . In the shell-model, the (e+PP)^^ decay modes of -^^Ca 
[25] . ^'^Ar, ^■^Fe, ^*Ni [26] and ^^Mo |27] isotopes have been investigated. 

The QRPA has emerged as the most successful model in explaining the observed 
quenching of NTMEs by incorporating the particle-particle part of the effective 
nucleon-nucleon interaction in the proton-neutron channel [551 [15] and the observed 
half-lives of several {PP)2i, decay emitters were reproduced successfully [7]. In 
the QRPA model, Staudt et al have evaluated the NTMEs for (/3+/3+)2,, mode 
[30] , The {e~^PP)2^ decay modes were studied in the QRPA and its extensions 
[13l[l4l[3Tl[32l[33l[i4l[35]. In spite of the success of the QRPA in the study of PP 
decay, there is a need to include the deformation degrees of freedom in its formalism. 
The deformed QRPA model has been developed for studying PP decay of spherical as 
well as deformed nuclei. However, these studies have been mainly restricted to P~ P~ 
decay so far [351 (SZl [35] ■ Besides the QRPA, the other models employed to study the 
{e+PP)2^ decay modes are the SU(4)^^ SSDH flT and pseudo-SU(3) [H]. 

In fact, the subtle interplay of pairing and quadrupolar correlations present in 
the effective two-body interaction decides the complex structure of nuclei. In addition 
to the pairing interaction, which plays an important role in all even Z-even N PP 
emitters, the crucial role of deformation degrees of freedom in the structure of ^^^Mo 
and ^^''Nd isotopes has been already reported [42l [43]. In the Projected Hartree- 
Fock-Bogoliubov (PHFB) model, the two crucial components of effective two body 
interaction, namely pairing and quadrupolar correlations are incorporated on equal 
footing and the rotational symmetry is restored by projection technique providing 
wave functions with good angular momentum for the parent and daughter nuclei 
involved in the PP decay. However, the PHFB model is unable to provide information 
about the structure of intermediate odd Z-oAA N nuclei in its present version and 
hence, on the single-/? decay rates and the distribution of GT strength. In spite of this 
limitation, the PHFB model, in conjunction with pairing plus quadrupole-quadrupole 
(PQQ) interaction [44] , has been successfully applied to study the 0+ 0+ transition 
of iP^P^)2^ modes |45l [Ml [47l [H] , where it was possible to describe the lowest excited 
states of the parent and daughter nuclei along with their electromagnetic transition 
strengths, as well as to reproduce their measured P~P~ decay rates [45l l48]. 

In the PHFB model, the existence of an inverse correlation between the 
quadrupole deformations and the magntudes of NTMEs M2,, has been shown [45l l47l 
[48] . In addition, it has been observed that the NTMEs are usually large in the absence 
of quadrupolar correlations. With the inclusion of the quadrupolar correlations, the 
NTMEs are almost constant for small admixture of the QQ interaction and suppressed 
substantially in realistic situation. For similar deformations of parent and daughter 
nuclei, the NTMEs have well defined maximum [13]. In the Interacting Shell Model 
(ISM) [50l[5T], similar observations have been also reported. Presently, we aim to study 
the e^PP decay of ^^^Dy isotope. The deformation parameters P2 for parent ^^^Dy and 
daughter i^^Gd nuclei are 0.2929 ±0.0016 and 0.3376 ±0.0018, respectively, which will 
provide a typical case study of the deformation effect. The present paper is organised 
as follows. We briefly present in section 2 the required theoretical formalism. In 
section 3, we calculate half-lives for the 0+ 0+ transition of {sP^)2^, and {ee)2^ 
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modes for ^^^Dy isotope together with various spectroscopic properties, specificahy, 
yrast spectra, reduced 5(^^2:0''' 2+) transition probabihties, quadrupole moments 
Q{2~^) and gyromagnetic factors 5(2"'") of ^^^Dy and ^^^Gd nuclei. The expressions 
to calculate these spectroscopic properties in the PHFB model are given in |52j . 
Moreover, we study the effect of deformation on NTME M21, of (e"*"/3/3)2j, decay modes 
vis-a-vis the changing strength of the QQ interaction. Finally, the concluding remarks 
are presented in section 4. 

2. Theoretical framework 

The inverse half-lives for the 0"'" 0^ transition of {e'^ PP)^^ decay modes are given 

by 



-1 



Tm(3)\ =G2Afi)\M2uf (1) 

where /3 denotes f3~^) 21,/ i^f^^) 21^/(^^)21' ™ode. The phase space factors G2u{(i) 
have been calculated to good accuracy ^ and the nuclear model dependent NTME 
M2V is expressed as 

To evaluate equation ([2]), one has to explicitly sum over all states of the intermediate 
odd Z-odd N nuclei. However, it is not possible to study the structure of intermediate 
odd-odd nuclei in the present version of the PHFB model. Alternatively, the 
summation over the intermediate states can be carried out by using the summation 
method [53j , when the GT operator commutes with the effective two-body interaction 
[Ml [551 mj. We use the latter procedure to sum over the intermediate states for 
evaluating the NTME [45l|471|48]. 

Presently, the nuclear wave functions are generated in the HFB framework by 
using a Hamiltonian consisting of an effective two-body interaction with pairing and 
quadrupolar correlations [44] . Explicitly, the Hamiltonian is written as 

H = H,p + V{P) + (:qqV(QQ) (3) 

where -ffsp, V{P) and V{QQ) denote the single particle Hamiltonian, the pairing and 
QQ part of the effective two-body interaction. Further, C^qq denotes the strength of 
QQ part of the effective two-body interaction. The purpose of introducing it is to 
study the role of deformation by varying the strength parameter C,qq. The final results 
are obtained by setting the C,qq = 1. 

With the assumption that the two-body part of this Hamiltonian almost 
commutes with the GT operator as the violation of isospin symmetry is presumably 
small, the NTME for the 0+ O"'" transition of (e+/3/J)2^ decay in the PHFB 
model in conjunction with the summation method is given by [471 148| 

^^^^ y (<^^°ll(^-^T-r-),J|vl/i(r") 



J.=OJf=0 1-1/2 
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where tt and v refer to protons and neutrons in the parent and daughter nuclei involved 
in the (e+/3/3)2^ decay modes. Further, 
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The expressions to calculate n"^ , n(^z,N),{z-2,N+2){(^)i fz-2.N+2 and Fz^n{0) are given 
in glllT]. 

In the present scheme, the difference in proton and neutron single particle energies 
with the same quantum numbers corresponding to the energy of the Isobaric Analog 
States is well described by the difference in Coulomb energy A^. Further, the spin- 
orbit energy splitting is added whenever required. The explicit inclusion of the spin- 
orbit splitting in the energy denominator of equation ([5]), implies that it cannot be 
factorized out of the sum in equation In the present context, it is noteworthy that 
the use of the summation method is richer than the closure approximation, because 
each proton-neutron excitation is weighted depending on its spin-flip or non-spin-flip 
character, and in addition, employing the summation method in conjunction with the 
PHFB formalism goes beyond what was done in previous applications of the pseudo- 
SU(3) model [51 [531 [H]. 



3. Results and discussions 

In the HFB framework, the nuclear wave functions for ^^^Dy and ^'"^^Gd isotopes are 
generated using the same model space and single particle energies (SPE's), but for 
the SPE's of O/111/2, 1/7/2 and 0/19/2 orbits, which are 4.6 MeV, 11.0 MeV and 11.6 
MeV, respectively, as in our earher calculation on (e+/3/3)2y decay modes [48]. The 
strengths of the pairing interaction are fixed as Gp = Gn = 30/A MeV. The strengths 
of the like particle components of the QQ interaction are taken as Xpp = Xnn = 0.0105 
MeV b~^, where b is oscillator parameter. For a given model space, SPE's, Gp, G„ 
and Xpp, we fix the strength of proton- neutron {pn) component of the QQ interaction 
Xpn by reproducing the excitation energy E2+ of the 2+ state. The adopted values of 
Xpn (in MeV b~^) for ^^^Dy and ^^^Gd isotopes in the present calculation are 0.01817 
and 0.02989, respectively. All these input parameters are kept fixed throughout the 
subsequent calculations. 

We present the calculated as well as experimentally observed results for yrast 
spectra [57], reduced _B(£'2:0+ 2+) transition probabilities [5^, static quadrupole 
moments Q(2+) and gyromagnetic factors g(2+) [59] of ^^^Dy and ^^^Gd nuclei 
in table 1. We tabulate only some experimentally observed representative values 
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in the case of reduced B{E2:0^ 2+) transition probabilities. The calculated 
B{E2:0~^ — > 2"*") are in excellent agreement with the observed ones. No experimental 
Q(2+) result is available for ^^^Dy isotope. The agreement between the calculated 
and experimental Q{2'^) for ^^^Gd nucleus is also quite good. However, the calculated 
gyromagnetic factors 5(2+) for ^^^Dy and ^^^Gd isotopes are off from the experimental 
data. 

The phase space factors of {sP^)2^ and {se)2^ modes for ^^^Dy are calculated 
following the notations of Doi et al 4 in the approximation Ci = 1.0, C2 — 0.0, 
C3 = 0.0 and i?i^i(£) = i?+i(e) + R-i{e) ~ 1.0. The calculated phase space factors 
are 4.723 x lO^^a y-i and G2U (ee) = 2.962 x lO-^o y-i. However, it is 

more justified to use the nuclear matter value of qa around 1.0 in heavy nuclei. Hence, 
the theoretical half-lives T'^'j'2 calculated for both gA — 1.0 and 1.261. In tabled 
we present the calculated NTME and half-lives Ty2 of {£P^)2u ^'^^ (^s)2i/ modes 
for ^^^Dy isotope along with other available theoretical results. No experimental 
result is available for the (e+/3/3)2^ decay modes of ^^^Dy isotope. Theoretically, the 
0+ 0+ transition of (ee)2y mode for ^^^Dy has been investigated only in the pseudo- 
SU(3) model gT]. The calculated NTME Ma^ in the PHFB model is 0.0138, which is 
smaller by a factor of 4.4 approximately than that due to pseudo-SU(3) model [41]. 
In the PHFB model, the calculated half-lives Ty2 '^^ i^(^~^)2u ^^'^ {^^)2v modes for 
gA = (1.261-1.0) are (1.115-2.819) x lO^^ y and (1.778-4.496) x 10^3 y, respectively. 

The multipolar correlations in general and quadrupolar correlations in particular 
present in the effective two-body interaction play a crucial role in the deformation of 
nuclei. Therefore, a natural choice to understand the role of deformation on the NTME 
M21, is to study its variation with respect to changing strength of the QQ interaction 
(qq. The results are presented in table [3] and figure 1(a). It is observed that the M21, 
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Figure 1. (a) Dependence of AI21, on the strength of QQ interaction {^qq. (b) 
NTMEs of (2/3+^^ and {ee)2^ modes for -"^^^Dy nucleus as a function of the 
difTcrencc in the deformation parameter A/32, "x" denotes the value of NTME 
for calculated A/32 at (^qq = 1. 

remains almost constant as the strength of (qq is changed from (qq ^ 0.0 — 0.4. As (qq 
is further increased up to 1.05, the NTME M21, starts decreasing except at (qq = 0.80, 
where a small increment in the size of M21, is noticed. For a further variation of Cqq 
from 1.05 to 1.50, the NTME M21, remains almost constant. In addition, a direct 
proportionality between the quadrupole moment Q2 and deformation parameter P2 
exists. In our earlier works [47l|48], it has been shown that the NTMEs A/2j/ are usually 
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large for (qq — 0.0 i.e. when both the parent and daughter nuclei are spherical. With 
the increase of (qq, the NTMEs remain almost constant and then decrease around 
the physical value Cqq = 1-0 establishing an inverse correlation between M21, and 
deformation parameter P2 ■ Presently, a similar inverse correlation between the NTME 
M21, and Q2 as well as P2 also exists. The effect of deformation on M21, is quantified 
by defining a quantity D2v as the ratio of M2U at zero deformation {(^qq = 0) and full 
deformation {(qq = 1). The ratio D21, is 13.64 for ^^^Dy nuclei. 

To investigate the observed suppression of the NTMEs as observed in the (3~ P~ 
decay with respect to the spherical case when the parent and daughter nuclei have 
different deformations [37l [49l [51], we present in figure 1(b), the NTME as 
a function of the difference in the deformation parameters A/32 = /32{parent) — 
P2{daughter) of the parent and daughter nuclei. The NTME is calculated by keeping 
the deformation of parent nucleus fixed at Cqq — 1 and the deformation of daughter 
nucleus is varied by taking (qq = 0.0 to 1.5. It is noticed from the figure 1(b) that 
the NTME is maximum when the absolute of the difference in deformation |A/32| is 
minimum and when the |A/32| increases the NTME M21, is strongly reduced. Thus, 
it is clear from this behaviour of M21, with respect to changing (qq, that the NTME 
tends to be large in the absence of quadrupolar correlations i.e. for a pair of spherical 
e+/3/9 parent and daughter /3/3 emitters. In figure [21 we present the 3-dimensional 
view of variation in NTME with respect to (qq — 0.0 — 1.5 for parent and daughter 
nuclei independently. The plateau corresponds to the maximum value of M21, for small 




1.5 



Figure 2. Dependence of Af2i/ on the independent variation of <^qq for parent 
and daugliter nuclei. 

admixture of quadrupolar correlations corresponding to I32{parent) — 0.0 — 0.086 and 
(32{daughter) — 0.0 — 0.163. The next peak is due to the increase in the magnitude of 
M21, around (qq — 0.8 and finally at C,qq — 1.0, the small physical M21, is obtained. 
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4. Conclusions 

To conclude, we have tested the rchabihty of HFB intrinsic wave functions for 
^^^Dy and ^^^Gd nuclei by calculating the spectroscopic properties, namely the 
yrast spectra, reduced B{E2:Q^ 2+) transition probabilities, static quadnipolc 
moments Q(2+) and 51-factors (?(2+) of these isotopes and comparing them with 
the available experimental data. An overall agreement between the calculated and 
observed spectroscopic properties suggests that the PHFB wave functions generated 
by fixing Xpn to reproduce the E2+ are quite reliable. Subsequently, we employ the 
same PHFB wave functions to calculate NTME as well as half-lives T^j2i^/3~^) 
and T^^2{£^) of '^^^Dy isotope for the 0+ 0+ transition. We also examine the effect 
of deformation on NTME M2V by varying the strength of QQ part of the effective two- 
body interaction. It is noticed that the M^v is the largest in the absence of quarupolar 
correlations. Moreover, it is reduced by a factor of 13.6 due to the quadrupolar 
correlations, which may be taken as a conservative estimate of the deformation effect 
in view of the approximations inherent in the present calculation. Further, the NTME 
has a well defined maximum when the deformations of parent and daughter nuclei are 
similar. 

This work was partially supported by DST, India vide sanction No. SR/S2/HEP- 
13/2006, by Conacyt-Mexico and DGAPA-UNAM. 
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Table 1. Excitation energies (in MeV) of = 2+ , 4+ and 6+ yrast states, 
reduced i3(i?2:0+ — > 2+) transition probabilities in b^, static quadrupole 
moments Q(2+) in e b, g-factors g{2'^) in nuclear magneton for ^^^Dy and ^^^Gd 
isotopes. Here B{E2) and Q(2+) are calculated for effective charge Cp =l+esff 
and en = £eff with e^ff = 0.60 and g(2+) has been calculated for g'j' = 1.0, 



= 0.0 and gj = 


g,- = 0.50. 














it-'^Gd 




Theory 


Experiment 


Theory 


Experiment 




0.1379 


0.1378" 


0.0886 


0.0889" 


^4+ 


0.4572 


0.4041" 


0.2936 


0.2882" 




0.9528 


0.7703" 


0.6119 


0.5847" 


B{E2) 


3.888 


3.710±0.040*' 


4.052 


4.64±0.05*'' 






3.46±0.33' 




3.6±0.5' 






3.79±0.30' 




4.5±0.5'' 


Q(2+) 


-1.786 




-1.823 


-1.93±0.04'= 










-1.96±0.04'= 


9(2+) 


0.550 


0.39±0.04'= 


0.602 


0.387±0.004'= 



"157J; "^f^ 



'denotes the adopted value. 



Table 2. Theoretically calculated and corresponding T^f (0+ 0+) of 

^£/3+^ and {e£)2„ modes for ^^^Dy nucleus. The numbers corresponding to (a) 
and (b) are calculated for gj\ = 1.261 and 1.0 respectively. * denotes the present 



work. 



Decay 
Mode 


Ref. 


Model 


\M2u\ 






eP+ 




PHFB 


0.0138 


(a) 


1.115x10^6 










(b) 


2.819xl02« 


ee 


* 


PHFB 


0.0138 


(a) 


1.778x10^3 










(b) 


4.496x10^3 




m 


p-SU(3) 


0.061 


(a) 


9.073x1021 










(b) 


2.294x10^2 



Table 3. Effect of the variation in C,qq on (Qg) , 02 and for ^^^Dy isotope. 



CI39 


ib 


•^Dy 


15« 


Gd 


A/2. 








{Ql) 


/32 




0.00 


0.000 


0.000 


0.000 


0.000 


0.191 


0.20 


0.093 


0.073 


0.071 


0.071 


0.187 


0.40 


0.345 


0.083 


1.672 


0.082 


0.187 


0.60 


28.33 


0.102 


65.81 


0.220 


0.015 


0.80 


65.98 


0.218 


80.82 


0.285 


0.029 


0.90 


74.38 


0.250 


86.39 


0.305 


0.022 


0.95 


82.82 


0.276 


88.72 


0.311 


0.018 


1.00 


91.60 


0.300 


90.78 


0.316 


0.014 


1.05 


94.27 


0.308 


92.10 


0.319 


0.003 


1.10 


94.88 


0.310 


92.58 


0.321 


0.002 


1.20 


95.42 


0.312 


93.25 


0.323 


0.003 


1.30 


96.08 


0.315 


93.66 


0.324 


0.003 


1.40 


96.62 


0.317 


94.00 


0.326 


0.002 


1.50 


97.01 


0.318 


94.32 


0.327 


0.002 



